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ABSTRACT 

Infection with high risk Human papillomavirus (HR-HPV) is necessary but not sufficient to cause cervical 
carcinoma. This study explored whether multiple HR-HPV or coinfection with Epstein-Barr virus (EBV) 
influence the integration status of HPV16 genome. The presence and typing of HPV in a series of 125 
cervical specimens were assessed by polymerase chain reaction (PCR) using the specific primers for the 
HPV LI region. As for EBV infection, the viral EBNA1 gene was used for its detection through PCR 
amplification. Disruption of the HPV E2 gene was assessed by amplification of the entire E2 gene with 
single set of primers, while E2 transcripts were evaluated by a reverse transcription PCR method (RT-PCR). 
The overall prevalence of HPVDNA was of 81.8% in cervical cancers versus 26.9% in benign lesions. In 
HPV positive cases, HPV 1 6 and HPV 1 8 were the most prevalent types, followed by HPV types 33,31. EBV 
EBNA1 prevalence was statistically more frequent in cervical carcinomas than in benign lesions (29.5%, vs 
9.6%; P=0.01). No viral infection was detected in healthy control women. The uninterrupted E2 gene was 
correlated with the presence of E2 transcripts originating from the HPV episomal forms. It was observed 
that integration was more common in HPV 18 and EBV coinfection. The presence of EBV caused a five-fold 
[OR= 5; CI= 1.15-21.8; P = 0.04] increase in the risk of HPV16 genome integration in the host genome. 
This study indicates that EBV infection is acting as a cofactor for induction of cervical cancer by favoring 
HPVDNA integration. 
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INTRODUCTION 

Human papillomavirus (HPV) infection is a prerequisite 
for the development of cervical cancer (6). This association is 



strong, consistent, and specific to a limited number of viral 
types (9, 26), in which HPV 16 and HPV 18 are the most 
prevalent (23, 14, 48). Geographic differences in the 
frequencies of other high-risk HPV types, including HPV 31, 
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33, and 45, have been reported to exist (10, 26). Nearly all 
women will become infected with HPV during their lifetime 
but only a minority of these infections will progress to invasive 
cancers. Therefore, there is a continuing need for the 
identification of viral and host factors that modulate the risk of 
disease progression (11). Possible agents included Herpes 
viruses among which more recently Epstein-Barr Virus (EBV), 
have been suspected to participate in cervical carcinogenesis 
(15, 33, 31, 13). For example, latent infection with EBV has 
been shown to act as a carcinogenic cofactor in several 
epithelial cell malignancies (38). The viral EBNA1 regulator of 
EBV transcription and replication is typically expressed in 
most forms of EBV-positive malignancies (24). On the other 
hand, integration of HPV DNA into the human genome is 
considered an essential step in the progression of HPV- 
associated cervical infection to invasive cervical carcinoma 
(18, 41). Upon integration the viral E2 gene, encoding a 
transcriptional repressor, is often disrupted with loss of its 
repressor function thus leading to up regulation of the 
transcription of the E6 and E7 oncogenes (30, 12). High 
expression of the viral oncogenes leads to disturbance of cell 
cycle controls and genomic instability that increase the risk of 
malignant transformation (25, 29). The main goal of our study 
was to explore frequency distribution of HR HPV genotypes 
and EBV infection in cervical cancer and premalignant lesions 
in order to correlate them with HPV 16 physical state. 

MATERIALS AND METHODS 

Sample Collection and Study Population 

This study enrolled 44 patients with cervical carcinoma 
recruited from the Radio-oncology department of Salah Azaiez 
Institute (Tunis, Tunisia). The mean patient age was 57.2 years 
(range: 38-76 years). Tumor specimens have been collected by 
punch biopsy, which consisted of 39 squamous cell carcinomas 
and 5 adenocarcinomas. Stages of the disease were coded 
according to International Federation of Gynecology and 



Obstetrics (FIGO). None malignant cervical biopsies were 
obtained from 52 women attending the Obstetrics and 
Gynecology Department of La Rabta Hospital (Tunis, Tunisia) 
for reasons other than cervical cancer. Also twenty nine 
Healthy controls women, counseling for routine reproductive 
healthcare, were used as reference. The range of ages varied 
from 27 to 56 years old, with an average of 41 years. Controls 
had no past or current history of cervical abnormalities. Biopsy 
specimens were transferred into cryotubes, suspended in 1 ml 
PBS (phosphate-buffered saline, pH 7.4) and stored frozen at - 
20°C until analysis. 

DNA and RNA isolation 

DNA extraction was performed using the QIAamp system, 
QIAamp DNA Mini Kit (Qiagen, Hilden, Germany), according 
to standard procedures as described in the kit, with detergent 
lysis buffer and proteinase K digestion. Total DNA was eluted 
in 100 ul of elution buffer and stored frozen at -20°C. RNA 
from frozen tissues was isolated using RNeasy Mini Kit 
(Qiagen, Hilden, Germany) as recommended by the supplier 
and dissolved in 30(xl RNase free water. 

The quality and quantity of the target nucleic acids were 
assessed by agarose gel electrophoresis and optical density 
measurement (absorbance at 260 nm). 

HR-HPV typing and EBV detection 

Purified genomic DNA was amplified using PCR 
technique for the HR HPVL1 and EBV EBNA1 genes 
(synthesized by GeneCust Europe, Dudelange, Luxembourg) 
using specific primers, as well as for the internal reference P- 
globin gene, to determine whether the isolated DNA was 
suitable for amplification. Thermal cycles were performed in a 
programmable PCR thermal cycler (Applied Biosystem, Foster 
City, USA), as previously described (46, 21). Each PCR 
experiment was performed with positive (HPV plasmids, 
cervical carcinoma previously proved to contain EBV) and 
negative (Blank water) controls. Samples tested positive for the 
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P-globin DNA control, confirming that none of the samples 
was inhibitory for PCR or lacked DNA. The Primer sequences, 



Primers for RT PCR assay 

HPV 1 6 E2 sense 5 ' -GAACTGC AACTAACGTTAGA-3 ' 

HPV16E2 anti- sense 5'-TCCATCAAACTGCACTTCCA-3' 
(3Actin sense 5 ' -AGCC ATGT ACGTTGCTATCC-3 ' 

P Actin anti-sense 5 ' -TTGGCGTACAGGTCTTTGC-3 ' 

Assessment of HPV16 E2 DNA and RNA integrity 

Adapting an approach previously described (4), the 
integrity of E2 gene was assessed using a single set of 
sequence-specific primers that were designed to amplify Open 
Reading Frame (ORF) fragments spanning the full length of 
HPV 16 E2 gene (Table 1). 

All positive samples for HPV 16 LI expression were 
subjected to PCR amplification of the E2 ORF region. The 
presence of E2 amplicon in agarose gel confirmed an intact E2 
gene, while lack of specific band indicated disruption. HPV 16 
positive samples were further subjected to amplification of E2 
transcripts (28). To control RNA integrity, RT-PCR reactions 
using P-actin specific primers were performed (34). One 
microgram total RNA was reverse-transcribed using the one 



amplicon sizes and cycling programs are illustrated in Table 1 . 



PCR program 

\°C l',58°C l',72°C l';X40 

t°C 1', 66°C l',72°C 1' ;X40 
i°C 1', 55°C l',72°C 1' ;X35 
i°C 1', 55°C l',72°C 1' ;X35 
i°C 1', 55°C l',72°C 1' ;X35 
5°C 1', 55°C l',72°C 1' ;X40 
CI', 58°C 40s,72°C 40s ; X 35 
C 1', 50°C2', 72°C 1'30;X35 

149 94°C40s,60°C l',72°C 1' ;X40 

500 94°C 30s, 50°C 30s, 72°C 1' ; X 30 

step RT PCR Kit (Qiagen, Hilden, Germany). The amplified 
products were electrophoresed on a 1.5% agarose gel, marked 
by a 100-pb DNA ladder (Gene Ruler™, Fermentas, Foster 
City, USA), stained with ethiduim bromide, and visualized 
under UV light. Images were recorded by using the Photo 
documentation system (Doc-Print II, VILBER LOURMAT, 
Marne-la-Vallee, France). 

Several precautions were taken to prevent false positive results. 
Different steps, such as sample preparation, nucleic acid 
extraction, amplification and post PCR, were performed in 
strictly separated rooms. 

Statistics 

Data were analyzed using SPSS, version 18.0 for windows 



Table 1. Polymerase Chain Reaction and Reverse Transcription primers, product length and programs. 



Primers for PCR assay 



Designation 



Product length (bp) 



HPV 16 LI sense 5'- 

HPV16L1 anti-sense 5'- 

HPV18L1 sense 5'- 

HPV18L1 anti-sense 5'- 

HPV31 LI sense 5'- 

HPV31 LI anti-sense 5'- 

HPV33L1 sense 5'- 

HPV33L1 anti-sense 5'- 

HPV45L1 sense 5'- 

HPV45L1 anti-sense 5'- 

PGlobin sense 5'- 

(3Globin anti-sense 5'- 

EBV EBNA1 sense 5'- 

EB V EBN A 1 anti-sense 5 ' • 

HPV 16 E2 sense 5'- 

HPV16E2 anti-sense 5'- 



GCAAGCAACAGTTACTGCGACGT-3 ' 
GCAACAAGACATACATCGACCGG-3 ' 
AAGGATGCTGCACCGGCTGA-3 ' 
CACGCACACGCTTGGCAGGT-3 ' 
CGTCCMARRGGA WACTGATC-3 ' 
TGTTTGTGCTGCAATTGCAAACAGTGATAC-3 ' 
CGTCCMARRGGA WACTGATC-3 ' 
TTTATGCACACAAGTAACTAGTGACAGTAC-3 ' 
CGTCCMARRGGA WACTGATC-3 ' 
ACACAAAATCCTGTGCCAAGTACATATGAC-3 ' 
CAACTTCATCCACGTTCACC-3 ' 
GAAGAGCCAAGGACAGGTAC-3 ' 
TGATAACCATGGACGAGGAC-3 ' 
CTTCAAGTTGCATTGGCTGC-3 ' 
ATGAAAATGATAGTACAGAC-3 ' 
CCAGTAGACACTGTAATAG-3 ' 



301 
217 
254 
111 
402 
268 
138 
1026 



94 
94 
95 
95 
95 
95 
94° 
95° 
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Vista. The associations between the variables were assessed 
using Fisher's exact test. Odds Ratios (OR) and corresponding 
95% Confidence Intervals (CI) were estimated by logistic 
regression models. A stepwise logistic regression analysis was 
used to assess the simultaneous effect of more than one 
variable and to identify for possible confounding factors. In all 
analyses, probability values less than 0.05 were regarded as 
significant. 

RESULTS 

Frequencies of HR-HPV DNA and EBV in cervical samples 

By using the primers directed against the highly conserved 
LI gene, we detected 4 different HPV types in our collection of 
samples. Overall 36 out of 44 cervical cancer samples (81.8%) 
were HPV positive. HPV type 16, the most frequent type of 
oncogenic HPV, was found in 24 cases (54.5%) while HPV18 
ranked second with a rate of 15.9 % (7/44). Less frequent types 
including HPV 33, and HPV31 were also found in 6.8 %, 4.5% 
of the samples respectively while HPV45 was not detected. 



Double infection was found in 5 cases of HPV 16 and 18 in 3 
cases and HPV 16 and 33 in 2 cases. 

From the benign cervical lesions, 21.1% were positive for 
HPV 16 (11/52), while 5.7% were positive for HPV 18 (3/52). 
One sample containing both HPV type 16 and 18 was found. 
Among the control group, no viral infection was detected. The 
results of HPV typing are summarized in Table 2. 

Our results indicate a significant increase in detection of 
HR-HPV in cancers compared to non malignant samples (P < 
0.01). 

The detection rate of EBV EBNA1 in tissue samples from 
cervical cancer patients was more frequent than that in non 
malignant samples (29.5% vs, 9.6 %, P = 0.01). In addition, 
EBV EBNA1 was not found in any of the control samples. 
According to the Odd Ratio (OR) [3.9 (1.28-12.16)], women 
infected with EBV have a three-fold likelihood to develop 
cervical carcinoma. 

A synergistic action of HPV 16 and EBV may also occur 
in some cervical carcinomas (22.7%), containing both EBV 
and HPV16 DNA versus only 7.6% in the benign samples. 



Table 2. Prevalence of HR-HPV types and EBV in cervical carcinomas and controls 





Controls 


Non malignant 


Cervical Cancer 


P values 


(N= 29) 


(N=52) 


(N=44) 


HPV types 


0 


14 (26.9) 


36(81.8) 




16 


0 


11 (21.1) 


24 (54.5) 




18 


0 


3 (5.7) 


7 (15.9) 




31 


0 


0 


2 (4.5) 


P<0.01* 


33 


0 


0 


3 (6.8) 




45 


0 


0 


0 




EBV 


0 


5 (9.6) 


13 (29.5) 


0.01* 


Multiple infections 


0 


5 (9.6) 


15(34) 




HPV16 & HPV18 


0 


1 (1.9) 


3 (6.8) 




HPV 16 & HPV33 


0 


0 


2 (4.5) 


0.004* 


EBV & HPV 16 


0 


4 (7.6) 


10 (22.7) 





Correlation between HPV16 E2 physical state and episomal 
derived transcripts 

Integrated forms of HPV 16 DNA were predominant in 
cervical cancers, constituting 62.5% (15/24), whereas episomal 
forms were detected in 9 cases (37.5%). Among benign lesions, 
2 of the 11 HPV 16 positive cases (18.1%) showed E2 



disruption, while episomal forms were noted in 9 cases 
(81.8%). Thus, E2 disruption was significantly higher among 
cancers as expected [P= 0.02, OR (95% CI) = 7.5 (1.32-42.77)] 
Table 3. We then decided to verify transcriptional expression 
of E2 in HPV 16- associated lesions using RT-PCR. 

All samples were found positive for the P-actin RNA 
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control, which was included to avoid false-negative results due 
to degradation of RNA (not shown). HPV16 E2 mRNAs were 
expressed in the presence of the intact virus with episomal 
forms and E2 transcripts originate only from the episomal 
HPV16 DNA. In Figure 1A and Figure IB, 6 of the 18 
HPV16E2 positive cervical cases (lanes 1-6) are shown as 
examples. 

Presence of the intact HPV16 E2 DNA, flanking the 1026- 



bp segments of the ORF fragment, was finally confirmed by 
HPV16 E2 RNA using primers giving an amplified product of 
149 bp. In general amplification of E2 with RT-PCR appeared 
more sensitive than amplification of the DNA-PCR method for 
detection of intact HPV E2 gene, because E2 DNA signals 
from RT-PCR are usually stronger than those from DNA-PCR 
(Figure 1 A, Figure IB). 



Table 3. E2-gene status in human papillomavirus type-16 positive benign lesions and cancers 



Parameters 


Total 

N 


Physical status of HPV16 DNA 

(%) 


OR (95% CI) 


P Value 




Episomal Integrated 






Benign cervical lesions 
Cervical cancers 


11 
24 


9(81.8) 2(18.1) 
9 (37.5) 15 (62.5) 


7.5 (1.32-42.77) 


0.02* 



P: two-sided Fisher's exact test 





1026 bp 

• Figure 1. Detection of HPV 16 physical status. (A) PCR 

301 bp products of HPV 16 episomal forms from cervical cancer 
positive samples: lanes 1-6 (amplicon length 1026 bp); M, Size 
markers of 100 bp ladder; N, Negative control (no template); 
C, positive control (HPV16 plasmid). 

(B) Transcripts product for positive HPV 16 E2 confirming 
PCR results: lanes 1-6 (amplicon length 149 bp); lane 1', 
. 1 indicative of integrated form witch was used as control. 



Disruption of the E2 gene together with multiple HR HPV 
or EBV coinfections predict carcinogenic progression 

To analyze whether viral factors might have an influence 
on the integration status, we correlated disruption of the 
HPV 16 genome with multiple HR HPV and EBV coinfections. 

Such analysis revealed that HPV 18 coinfection was found 
more often associated with HPV 16 integration, while in 
contrast, HPV33 coinfection was only associated with HPV 16 



episomal forms (Figure 2). We also correlated the presence of 
EBV to the status of the HPV 16 genome in coinfected lesions. 
Interestingly, EBV was more frequently present in samples 
with integrated forms of HPV 16 than in the episomal forms. 
Based on these data, EBV infection was significantly 
associated with disruption of the HPV16 genome as revealed 
by OR (95% CI) and P values [OR= 5; CI= 1.15-21.8; P = 
0.04] Table 4. EBV could therefore be considered as a cofactor 
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of disease risk as determined by the OR based on logistic In conclusion of our analysis, EBV may play a 

regression models. The presence of EBV also may increases cooperative role in HPV16 carcinogenic progression by 
the probability of integrated HPV16 genome. favoring HPV16 DNA integration in the cellular genome. 




e 2. Frequency distribution of HR 
types and EBV correlated with 
6 DNA physical status. 



KJV1J WIS PIN BV 



Table 4. Distribution of EBV positive and negative cervical specimens correlated with HPV 16 DNA physical status 
HPV16 Physical Total EBV 

nrviornvMuu lo ai ^ OR (95% CI) P Value 

status N Positive Negative 

T E fTi !5 *fSl3 5(1.15-21.8) O^H 
Integrated 17 10 (58.8) 7 (41.1) ^ _ 

P: two-sided Fisher's exact test 



DISCUSSION 

Infection with HR HPV is a necessary but not sufficient 
cause of cervical carcinoma, with additional viral and host 
genetic events required to drive cells to the malignant 
phenotype (47). In this report, we confirmed that HPV16 and 
18 are the two most frequently detected types among cervical 
carcinoma in Tunisia, and that they account for 70.4%, which 
is very similar to the proportion estimated in other world 
regions by recently published large Meta analyses (23, 14). To 
a lesser extent, HPV31 and HPV33 were detected in 11.3 % of 



all cases, whereas HPV45 was not found. Similar data have 
been described for HPV31, 33, 35, 45, 52 and 58, which can be 
grouped in several branches differing in geographic 
distribution, and in relative prevalence within different ethnical 
groups (7, 8). HPV 18 is the primary HPV type associated with 
cervical adenocarcinomas and its related type HPV45 (alpha 7 
species) shows a similar rate of association (42). Interestingly, 
our pathological data demonstrate larger predominance of 
squamous cell carcinoma over adenocarcinomas in our 
samples, Furthermore the presence of HPV in benign cervical 
lesions has been widely documented. In our study, type 
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specific HPV prevalence was similar to that found in 
aggregative analysis, where HPV 16/18 are the most frequent 
types observed (39, 40). 

Although the role of high risk HPV infection in cervical 
cancer is well documented, the role of multiple infection has 
been less studied. Multiple HPV infection has been found to be 
associated with persistence of HPV infection that is a 
prerequisite for neoplastic transformation (2, 17). That report 
detected multiple HPV infection in 11.3% of cervical cancer 
patients. In other reports, the rate of multiple HPV infection in 
invasive cervical cancers ranged between 0 % and 36 % (32). 
Thus, multiple HR HPV infections ought to be more relevant 
for HPV18. Women coinfected with HPV 16 and HPV 18 are at 
significantly higher risk of developing CIN3 and cancer when 
compared to women infected with other high-risk types (22). 

Interestingly, infection with HPV 33 was almost 
exclusively detected as part of multiple HPV infection. Only 
one patient had HPV type 33 as a single infection among the 
three cases detected. One explanation as to why HPV type 33 
was a common element of multiple infections may be that HPV 
33 usually occurs in combination with phylogenetically closely 
related types (HPV 16, 31, 35 and 58). On the other hand, 
HPV3 1 was considerably less common in cervical carcinoma, 
as seen in previous studies (20), suggesting that in Tunisian 
populations, it has only a relatively weak carcinogenic 
potential in comparison to other types. Only one non malignant 
lesion was shown to harbor HPV16/18 coinfection. It is likely 
that this specific benign lesion is at risk of carcinogenic 
progression. 

Otherwise, prevalence of HR HPV types alone does not 
give an accurate picture of the carcinogenic state of the lesion. 
In this regard, studies of the physical status of the HPV genome 
could help evaluating the risk of cervical carcinogenesis. Our 
findings confirmed earlier reports on association of disruption 
of the E2 gene with cervical cancer, wherein squamous cell 
carcinomas harbored 62.5% disruption in contrast to 18.1% in 
benign lesions. These data are in agreement with the report by 



Nagao et al., (27), who showed an increase of integrated forms 
of HPV 16 in parallel with the disease progression. The HPV 
E2 protein plays a crucial role in the viral life cycle as a 
transcription and a replication factor. E2 was also shown to be 
able to shuttle from the nucleus, to the cytoplasm where it 
activates caspase 8, thus inducing apoptosis (5, 36, 37). The 
HPV 16 E2 protein has also been shown to be expressed at 
relatively high levels in benign lesions compared to cancers 
(45). Consequently, the detection of HR-HPV E2 transcripts 
might serve as an additional biomarker of cervical diseases. 
Results from the current study also indicate that the 
transcriptional activity of HPV16 E2 is a more sensitive 
indicator of the genomic disruption than analyses of the 
genomic DNA. When cervical carcinomas coinfected by 
HPV 16 and HPV 18 were investigated we could show that 
HPV 18 coinfection was more likely linked to HPV 16 
integrated forms. It was previously established that HPV 16 is 
present exclusively in episomal forms in 30-70% of cervical 
cancer (3, 43), whereas HPV18 has been reported to be mainly 
integrated (3), which indicates that these HPV types might 
have different biological characteristics. Meanwhile, some data 
have indicated that patients infected with HPV 18 showed a 
rapid progression through precancerous stages (44). In contrast, 
HPV33 was found substantially more frequently associated 
with HPV16 episomal forms, in comparison with HPV18. It is 
in good agreement with the observation that cancers associated 
with HPV 16 and HPV 18 tend to be more aggressive, and that 
the detection of HPV types other than 16 and 18 is an 
independent predictor of better survival in patients with 
cervical carcinoma (19). Portraying the presence of intact E2 
genes in almost 37.5% of HPV 16 positive cervical cancer 
cases, have subsequently paved the way for new paradigms of 
cervical carcinogenesis. 

A putative role of cooperative factors such as interaction 
with co infecting viruses, possibly at different stages of 
cervical carcinogenesis is demonstrated here. Previously, EBV, 
an oncogenic herpes virus, has been suggested as a possible 
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viral cofactor (35). The fact that the EBV/C3d receptor was 
shown to be expressed in cervical epithelial samples, as well as 
the clear EBV linkage with epithelium originated tumors, made 
this hypothesis plausible (15). However, other authors did not 
reach to the same conclusion, confirming that EBV was acting 
in cervical tumor pathogenesis (33). The absence of EBV 
infection in controls and the decline of EBV peak in benign 
cervical lesions suggest that EBV infection occurs late in 
cervical oncogenesis. It may be that the EBV virus is acting as 
a cofactor of HPV, in induction of uterine cervix pathology; the 
suggestion was confirmed by studies of Szkaradkiewicz et al., 
(35), which report a possible sexual route for transmission of 
EBV. However, there is no convincing evidence to support the 
hypothesis that EBV infection plays a direct role in the 
pathogenesis of cervical neoplasia. Though, EBV-infected 
tumour infiltrating lymphocytes might contribute indirectly to 
cervical carcinogenesis via producing viral IL-10 expressed 
from the BamHI-C fragment rightward reading frame 1 gene 
(BCRF-1), which causes a reduction in local immunity leading 
to suppression of the response to HPV-transformed cells (1). 
This, together with the fact that EBV infection may contribute 
in the integration of HPV16 genome occurrence, suggests that, 
although they might not be involved in the oncogenic processes 
directly, they might play a role as cofactor. Accordingly, the 
viral EBV EBNA1 protein has been shown to be able to 
decrease apoptosis and DNA repair in EBV-associated lesions 
thus potentially contributing to malignant transformation of 
cotransfected EBV/HPV lesions (16). 

In summary, further knowledge of markers of HPV 
genomic integration has the potential to improve objective 
management of cervical neoplasia. The data presented here 
might also be beneficial to our understanding of the biology of 
human neoplasms associated with integration of other DNA 
tumor viruses. 

List of abbreviations: EBV: Epstein-Barr virus; RT-PCR: 
reverse transcription polymerase chain reaction; PBS: 
phosphate -buffered saline; HPV: human papillomavirus. 
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